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ABSTRACT 


The  propulsion  performance  of  a  submerged  body  of  revolution  was 
determined  when  the  axial  clearance  between  a  stern  propeller  and  the  body 
was  varied.  Theoretical  and  experimental  results  indicate  that  variations, 
within  the  range  of  investigation,  in  the  longitudinal  position  of  a  stern  pro¬ 
peller  do  not  critically  affect  the  propulsive  coefficient.  In  addition  to  these 
results,  calculation  of  the  propellers  and  the  thrust  deduction  is  discussed, 
and  experimental  results  obtained  from  numerous  wake  traverses  are  presented. 


INTRODUCTION 

As  part  of  the  Bureau  of  Ships  Fundamental  Hydromechanics  Research  Program,  the 
David  Taylor  Model  Basin  is  studying  the  problem  of  propulsion  interaction.  Propulsion 
interaction  effects  for  a  single  propeller  and  for  counterrotating  propellers  on  deeply  sub¬ 
merged  bodies  of  revolution  have  been  the  subject  of  previous  investigations.1'2'3*4 
Weinblum5  in  his  paper  on  thrust  deduction  has  concisely  set  forth  the  physical  nature  of 
the  propulsion  interaction  problem.  In  a  general  manner  he  described  some  of  the  theoretical 
methods  used  in  calculation  of  thrust  deduction  and  indicated  the  state  of  recent  scientific 
knowledge.  Breslin  showed,  in  Reference  2,  the  dependence  of  the  thrust  deduction  upon  the 
slope  of  the  sectional  area  curve  at  the  stern  of  a  body  of  revolution.  He  also  found  the 
variation  of  thrust  deduction  (for  a  given  size  propeller)  with  axial  clearance  for  an  ellipse 
of  revolution  having  the  same  fineness  ratio  as  a  hypothetical  arrangement  of  the  airship 
AKRON.  Korvin-Kroukovsky 1  calculated  the  wake  fraction  and  thrust  deduction  for  AKRON 
when  a  hypothetical  stern  propeller  installation  is  specified.  Historically,  a  great  deal  more 
experimental  work  has  been  done  on  the  propulsion  interaction  problem  as  associated  with 
surface  ships  than  has  been  done  on  submerged  bodies.  However,  considerable  theoretical 
work  has  been  accomplished  which  deals  with  the  various  aspects  of  the  interaction  problem 
for  the  case  of  axisymmetric  flow.  Emphasis  on  theoretical  work  in  this  direction  is  because 
the  axisymmetric  flow  case  is  more  amenable  to  mathematical  treatment. 

Each  of  the  aforementioned  papers  has  contributed  considerably  to  the  available 
information  on  the  subject  of  marine  propulsion;  however,  the  present  study  differs  in  several 
particulars  from  the  investigations  just  cited.  In  addition  to  the  determination  of  thrust 
deduction  and  wake,  the  total  propulsive  coefficient  is  considered  with  emphasis  on  obtain¬ 
ing  the  best  propeller  efficiency  for  a  given  set  of  operating  conditions.  This  is  extended 
to  give  the  performance  of  propellers  as  related  to  the  axial  spacing  between  propeller  and 
body.  A  Series  58  form  was  used  for  the  experiments  and  the  computed  example.  This 
report  includes  model  propulsion  test  results  and  presents  and  discusses  necessary  wake 
surveys,  thrust  deduction  determination,  open-water  data,  and  propeller  calculations. 


^•firwcei  are  listed  on  page  47. 


STATEMENT  OF  PROBLEM 


Propulsion  characteristics  of  a  submerged  body  are  determined  by  the  performance  of 
a  body  and  propeller  system.  The  body  affects  the  propeller  performance  through  its  wake 
and  the  body  resistance  is  affected  by  the  induced  velocity  field  of  the  propeller.  The  general 
problem  is  to  find  in  what  way  the  propeller  is  influenced  by  the  body,  and  conversely,  in 
what  way  the  body  is  influenced  by  the  propeller.  In  this  report  these  mutual  influences  are 
investigated  for  a  well-streamlined  body  of  revolution  when  equipped  with  a  stern  propeller 
located  at  different  axial  positions.  The  specific  problem  treated  is  to  determine  the  best 
efficiency,  as  a  function  of  axial  position,  of  a  stern  propeller  of  a  given  diameter,  and  then 
to  determine  the  optimum  axial  spacing  between  propeller  and  body  based  on  propulsive 
coefficient. 

The  approach  used  in  this  leport  to  reach  a  solution  to  this  interaction  problem  may 
be  outlined  as  follows:  The  velocity  and  pressure  distribution  over  the  surface  of  an  axi- 
symmetric  body  without  propeller,  in  potential  flow,  was  obtained  by  generating  the  pre¬ 
scribed  shape  with  a  surface  distribution  of  sources.  The  induced  velocity  field  ahead  of 
a  propeller  was  determined  by  assuming  that  the  propeller  may  be  represented  by  a  cylindri¬ 
cal  vortex  sheet.  The  thrust  deduction  force,  i.e.,  interaction  force,  was  obtained  by  suitable 
integration  of  the  pressure  distribution  over  the  afterbody  with  and  without  a  propeller.  The 
inflow  to  the  propeller  is  mostly  viscous  in  origin;  therefore,  the  wake  was  experimentally 
determined  at  various  positions  behind  a  Series  58  form  of  7.3  fineness  ratio.  An  optimum 
wake-adapted  propeller  was  designed  for  each  axial  position,  using  a  method  in  current  use 
at  the  Model  Basin.  To  investigate  the  effect  of  axial  position  on  propeller  performance, 
the  propeller  thrust  and  diameter  were  kept  constant.  Finally,  the  approach  as  outlined  per¬ 
mitted  the  propulsive  coefficient  to  be  computed.  Experimental  confirmation  was  obtained 
from  tests  of  one  calculated  propeller. 

Other  limitations  and  assumptions  involved  in  the  method  may  be  stated  as  follows: 

1.  It  is  assumed  that  the  propeller-induced  velocity  normal  to  the  body  can  be  neglected; 
i.e.,  the  boundary  conditions  on  the  body  are  not  met.  It  is  shown  that  the  accuracy  of  the 
calculated  thrust  deduction  is  not  seriously  impaired  by  this  neglect  for  the  well  streamlined 
body  considered. 

2.  It  is  assumed  that  the  fluid  is  incompressible  and  inviscid  except  that  the  total  wake 
to  the  propeller  is  considered;  i.e.,  the  problem  is  treated  as  a  potential  one  except  that  a 
viscous  correction  is  made  to  the  wake.  The  problem  must  be  treated  this  way  since  the 
effect  of  a  stern  propeller  on  a  body  of  revolution  is  mainly  potential  in  origin  whereas  the 
effect  of  the  body  on  the  propeller  is  mostly  viscous  in  origin. 
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ANALYSIS 


PRINCIPAL  CONSIDERATIONS 

To  make  the  theoretical  analysis  tractable,  only  a  submerged  well-streamlined  body  of 
revolution  was  considered.  This  choice  of  body  shape  does  not  infer  that  the  effects  of 
various  appendages  and  appendage  arrangements  are  not  of  practical  importance. 

From  a  propulsion  viewpoint,  the  performance  of  any  given  body  and  propeller  system 
can  be  analyzed  by  means  of  the  propulsive  coefficient  and  its  components.  The  relation 
is 


Vp  35  -p-  meh  *  ep  *  er 
1  -  t 

where  eh  is  the  hull  efficiency  — — — —  , 

1  -  Ir  o 

ep  is  the  propeller  open-water  efficiency,  and 
er  is  the  relative  rotative  efficiency. 

Consider  the  components  of  the  propulsive  coefficient.  First,  the  hull  efficiency  eh  is  the 
1  -  t 

ratio  - — —  ,  where  (1  -  ^0)  is  the  effective  wake  factor  which  is  experimentally  deter- 

1  “  "  o 

mined  and  (1  -  t)  is  the  thrust  deduction  factor.  It  is  necessary  to  obtain  wake  experimental¬ 
ly  inasmuch  as  the  theory  of  free  turbulent  shear  flow  is  not  developed  sufficiently  to  give 
reliable  results  in  the  near  field  behind  a  body.  The  flow  field  of  the  body  alone  is  regarded 
as  composed  of  two  superimposed  fields.  These  two  fields  are  the  body  velocity  V  and  the 
disturbance  velocity  W0V.  If  the  velocity  field  of  a  propeller,  which  is  operating  behind  a 
body,  is  superimposed  on  the  velocity  field  of  the  body,  the  combined  velocities  produce  a 
pressure  gradient  which  usually  reduces  the  pressure  along  the  afterbody.  This  pressure 
reduction  augments  the  resistance  of  the  body  so  that  only  part  (1  - 1)>  the  thrust  deduction 
factor,  of  the  propeller  thrust  is  available  to  overcome  the  resistance  experienced  by  the 
body  without  propeller.  The  thrust  deduction  problem  is  analogous  to  the  familiar  problem 
of  finding  the  forces  between  two  bodies  in  a  flow  field.  Unlike  the  inflow  into  the  propel¬ 
ler,  the  induced  velocity  field  ahead  of  the  propeller  and  the  body  velocity  without  a  propel¬ 
ler  can  be  computed  theoretically  with  sufficient  accuracy  for  the  present  investigation. 
Alternatively,  a  good  approximation  to  the  thrust  deduction  can  be  obtained  experimentally 
by  means  of  stock  propellers.  Experimentation  of  this  type  is  costly,  and  inasmuch  as  an 
adeqtfhte  theory  appears  to  exist  for  axisymmetric  bodies,  it  seems  preferable  to  seek  a 
theoretical  solution  with  experimental  confirmation  at  a  few  points* 


♦The  propul. ive  coefficient  should  not  be  used  to  compare  different  body  shape.. 
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Second,  consider  propeller  efficiency  behind  a  body.  It  is  apparent  that  the  calculation 
of  an  optimum  propeller  for  each  axial  position  behind  a  body  of  revolution  requires  the  use  of 
(1)  the  wake  which  is  obtained  experimentally  and  (2)  the  thrust  deduction  which  is  obtained 
theoretically.  The  circulation  theory  as  applied  to  wake-adapted  propellers  is  used  to  compute 
the  propeller  performance  at  each  axial  position.  Propeller  efficiency  for  the  behind  condition 
is  given  by  the  product  of  open-water  efficiency  ep  and  the  so-called  relative  rotative  efficiency 
ef .  If  a  wake  suney  is  available,  the  propeller  efficiency  for  the  behind  condition  ( ep  •  er)  may 
be  estimated  from  the  relation 


Next,  consider  propeller  performance  in  open  water.  Propeller  performance  character¬ 
istics  in  open  water  are  of  considerable  interest  particularly,  for  systematic  series  and  for 
ship  powering  estimates  where  the  concept  of  an  equivalent  speed  of  advance  is  introduced. 
An  effective  wake  fraction  W0  is  determined  by  this  procedure.  Making  use  of  an  effective 
wake  fraction,  propeller  open-water  efficiency  of  a  wake-adapted  propeller  can  be  approxi¬ 
mated  by  the  following  equation: 


where  the  propeller  thrust-load  and  power  coefficients  C7Si  and  Cp$i  are  based  on  ship 
speed  and  are  for  nonviscous  flow, 

T  is  propeller  thrust, 

Q  is  propeller  torque, 

Va  is  propeller  speed  of  advance, 

n  is  propeller  frequency  of  revolution, 

«  is  the  section  drag-lift  ratio, 
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fij  is  the  hydrodynamic  pitch  angle,  and 
R 

—  is  the  radius  fraction. 

*0 

Finally,  the  relative  rotative  efficiency  is  obtained  numerically  from  Equations  [1]  and  [la] 
as  the  ratio  of  the  propeller  efficiency  in  the  behind  condition  to  that  in  open  water. 

The  subsequent  sections  of  the  report  are  chiefly  concerned  with  the  details  of  reach¬ 
ing  a  solution  for  and  calculating  the  variation  of  the  propulsive  coefficient  with  axial  position 
of  a  propeller.  The  foregoing  approach  used  to  find  the  propulsive  components  is  summarized 
as  follows: 


Problem 


Method  of  Solution 


Determination  of  flow  due  to  the  body 

Determination  of  the  flow  ahead  of  a 
propeller 

Thrust  deduction 


Propeller  efficiency 


Experimental  wake  measurements  and  calcu¬ 
lation  of  pressure  distribution  on  body. 

Induced  velocities  due  to  a  cylindrical  vortex 
sheet. 

Body  and  propeller  flow  combined.  The  appropri¬ 
ate  pressures  on  the  hull  afterbody  then  inte¬ 
grated  to  yield  the  thrust  deduction  force. 

Circulation  theory  as  applied  to  wake-adapted 
propellers  using  a  design  method  developed  by 
Eckhardt  and  Morgan. 


DETERMINATION  OF  THRUST  DEDUCTION 

The  addition  of  the  propeller-induced  velocities  to  the  velocities  at  the  body  produces 
a  reduction  in  pressure  on  the  body  within  the  field  of  propeller  influence.  The  difference 

£iRt  «  T  -Rt 

where  T  is  the  propeller  thrust  and  R t  is  the  body  resistance  without  a  propeller,  is  variously 

A  Rt 

called  the  thrust  d3^uction  force  or  resistance  augmentation  and  the  ratio  t  =*  — —  is  known 

as  the  thrust  deduction  coefficient.  The  thrust  deduction  will  be  calculated  in  this  report  by 
pressure  integration.  The  formulation  of  the  solution  is  two-fold:  (1)  A  propeller  field  as 
represented  by  Korvin-Kroukovsky1  will  be  used  to  obtain  the  propeller-induced  velocities. 
The  component  of  these  velocities,  which  is  tangent  to  the  body,  is  then  combined  with  the 
potential  velocities  ut  at  the  body  surface  to  yield  the  total  fluid  velocity  V .  (2)  The  pres¬ 
sure  defect  A p  on  the  body  due  to  the  working  propeller  is  calculated  and  then  integrated  in 
the  z-direction  over  the  afterbody  to  yield  the  thrust  deduction. 

Turn  now  to  the  problem  of  combining  the  propeller  velocity  and  body  velocity  fields. 
Figure  1  diagrams  the  disturbing  velocities  due  to  the  propeller  and  shows  schematically  a 
typical  body  meridian  profile  and  the  propeller  disk  at  several  locations.  The  longitudinal 
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Figure  1  -  Schematic  Diagram  of  Body- Propeller  Arrangement  Showing 
th8  Magnitude  of  Pertinent  Velocities 

position  of  each  propeller  plane  is  a  distance  Lq  aft  of  a  so-called  reference  station  at 
L0  =  0.  For  the  remainder  of  this  report,  L0  is  expressed  as  a  fraction  of  the  maximum 
diameter  D  of  the  body.  The  magnitude  of  the  propeller-induced  velocity  u.  which  is  tangent 
to  the  body  is 

ui  m  vr  cos  (0  -  5)  -  vq  sin  (6  -  8) 
and  the  total  fluid  velocity  V  for  the  body  plus  propeller  is 

V  *  +  ut 

where  «f.  is  the  propeller-induced  velocity  tangent  to  the  body  surface, 
ut  is  the  potential  velocity  on  the  body  surface  without  propeller, 

V  is  the  total  fluid  velocity, 

vr  is  the  disturbing  velocity  directed  opposite  to  the  radius  vector  from  the  propeller, 

vq  is  the  disturbing  velocity  normal  to  the  spherical  radius  vector  (positive  when 
directed  toward  the  body  axis), 

6  is  the  vectorial  angle,  and 

8  is  the  angle  of  inclination  of  the  tangent  to  a  body  of  revolution. 

When  $  <  5,  vq  increases  V  and,  conversely,  when  6  >  5,  vq  decreases  V .  The  total  value  of 
V  for  propeller  and  body  must  be  greater  on  the  afterbody  than  the  value  of  U  for  the  body 
alone;  i.eM  the  increase  in  velocity  over  the  hull  afterbody  decreases  the  absolute  pressure. 

A  procedure  for  obtaining  the  total  fluid  velocity  ratio  V/V  is  formulated  in  detail  in  Appendix  C. 


e 


NUMERICAL  COMPUTATIONS  AND  RESULTS 


THRUST  DEDUCTION  FOR  CHOSEN  BODY 

The  pressure  distribution,  which  corresponds  to  a  source-sink  distribution  for  Model 
4620  without  a  propeller,  and  the  combined  pressure  distributions  of  this  body  plus  propeller 

are  shown  in  Figure  2.  Calculations  leading  to  the  distribution  of  the  pressure  field  A  — 
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in  Figure  2  are  arranged  in  Tables  4,  5,  and  6  of  Appendix  C  for  propeller  positions  —  of 

1/4,  1/2,  and  1,  respectively.  It  can  be  seen  from  Figure  2  that  the  effect  of  the  propeller 
field  is  significant  for  approximately  the  last  20  percent  of  the  body  length.  Propeller-induced 
velocities  normal  to  the  body  are  neglected;  thus,  the  boundary  condition  that  the  normal 
velocity  at  a  point  on  the  body  surface  be  zero  is  not  satisfied.  For  the  present  case,  the 
neglect  of  this  component  does  not  seriously  impair  the  accuracy  of  the  calculated  thrust 
deduction.  The  following  data  illustrate  the  order  of  magnitude  of  the  pertinent  quantities. 


Propeller-induced  velocity  vq  is  seen  to  be  essentially  zero  and  propeller-induced  velocity 
vr/V  sin  (0  -  3),  which  is  normal  to  the  tangent  line  at  a  point  on  the  body,  is  also  seen  to 
be  negligible.  When  sin  (0  -  8)  is  significant,  then  vr  is  small.  Of  course,  these  compen¬ 
sating  effects  are  a  result  of  the  streamlined  shape  of  the  body  under  consideration  and  may 
not  be  true  in  general.  The  integral  in  Equation  [4]  of  Appendix  C  can  be  computed  by 
Simpson’s  rule  using  the  slope  data  of  Table  3  (see  Appendix  C)  and  the  pressure  coeffi¬ 
cients  of  Figure  2.  Results  of  evaluating  t  as  a  function  of  LQ/D  are  shown  in  Figure  3. 

In  retrospect  the  wake  and  thrust  deduction  are  mutually  dependent  quantities.  As¬ 
suming  the  propeller  is  an  actuator  disk,  the  exact  functional  relation  between  t  and  W  in 
pure  streamline  flow  has  been  independently  derived  by  several  authors. 5»6»7  The  relation 
can  be  expressed  in  the  following  convenient  form: 

1-wp’1~  +  [2] 

The  potential  velocity  ratio  (1  -  Wp)  was  computed  from  Equation  [2]  by  making  use  of  the 
values  of  t  given  in  Figure  3  and  the  nonviscous  values  of  C T  which  were  computed  in 
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Figure  2a  —  Differences  between  Pressure  Coefficients 
for  Towed  and  Propelled  Conditions 


Figure  2b  —  Pressure  Distribution  with  and  without  a  Propeller 


Figure  2  -  Pressure  Field  on  Afterbody  of  Model  4620 
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Figure  3  -  Thrust-Deduction  Coefficient  versus  L0/D 


connection  with  the  propeller  design.  Results  are  plotted  in  Figure  21  of  Appendix  D.  These 
comparative  data  reveal  clearly  the  effect  of  nonuniform  wake  and  the  convergence  of  the 
curves  at  large  values  of  L0/D. 

GEOMETRY  AND  DESIGN  CONDITION  FOR  EACH  PROPELLER 

All  propellers  calculated  and  tested  are  five-bladed  designs  having  sections  comprised 
of  a  thickness  form  with  an  NACA  66  nose  and  a  parabolic  tail  with  an  a  =  0.8  mean  line. 

Each  propeller  was  calculated  by  the  circulation  or  vortex  theory  using  the  design  method  of 
Eckhardt  and  Morgan.8  The  approach  was  to  determine  the  optimum  propeller  from  the  stand¬ 
point  of  efficiency  for  each  longitudinal  position.  Although  cavitation  is  not  part  of  this  study, 
incipient  cavitation  curves  were  used  to  help  determine  the  blade-section  chord.  To  system¬ 
atically  investigate  the  effect  of  axial  spacing  between  body  and  propeller,  the  most  important 
remaining  parameters  should  remain  constant;  therefore,  to  satisfy  this  requirement  both  pro¬ 
peller  thrust  and  propeller  diameter  were  fixed.  Calculation  of  a  propeller  for  a  specified 
longitudinal  position  was  based  on  the  radial  wake  distribution  (see  Appendix  D)  and  the 
theoretically  obtained  thrust-deduction  coefficient  at  that  particular  transverse  station,  a 
constant  thrust  of  82,630  lb  for  a  200-ft  prototype  and  a  propeller- body  diameter  ratio  d/D  of 
0.5.  The  choice  of  0.5  for  the  value  of  d/D  is  based  on  TMB  classified  data  which  indicate 
that  d/D  «  0.5  gives  the  proper  propeller  diameter  in  the  present  case  for  the  behind  condition. 

A  significant  departure  from  what  has  been  done  by  most  investigators  in  dealing  with 
the  interaction  problem  is  the  emphasis  placed  on  a  suitable  propeller.  In  this  study,  an 
optimum  propeller  and  a  nonoptimum  propeller,  for  a  given  diameter,  are  determined  for  the 
flow  and  load  conditions  that  exist  at  each  location.  The  design  conditions  can  now  be 
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formulated,  based  on  the  values  of  (1-0  and  (1  -W  ). already  determined.  The  optimum  rpm- 
diameter  relationship  was  determined  from  Gawn's9  open-water  curves.  The  results  are: 


P/d 

J 

SI 

'^optimum 

1/8 

1.44 

1.16 

154.4 

1/4 

1.51 

1.20 

153.7 

1/2 

1.60 

1.27 

150.8 

1.0 

1.70 

1.36 

147.0 

Alternatively,  ep  and  P/d  were  obtained  for  a  range  of  rpm,  and  ep  is  plotted  in  Figure  4. 

For  comparative  purposes,  a  set  of  four  nonoptimum  propellers  which  would  operate  at  a 
higher  and  possibly  more  practical  rpm  was  calculated.  These  propellers  were  calculated 
for  200  rpm,  where  from  Figure  4  it  is  seen  that  a  drop  in  ep  of  about  6  or  7  points  has 
occurred.  The  loss  in  ep  is  much  more  rapid  in  the  range  of  N  >  200  and  demonstrates 
clearly  the  penalty  in  propeller  efficiency  for  operation  at  higher  rpm  for  the  chosen  di¬ 
ameter.  Altogether,  the  performance  of  nine  propellers  was  theoretically  determined,  one  was 
built  and  tested.  Table  1  summarizes  the  design  conditions  and  location  of  each  propeller. 

Calculations  for  moderately  loaded 


TABLE  1 

Summary  of  Propeller  Locations 
and  Design  Conditions 

T  =  82,630  —=  0.5  z  =  5 

D 


^0 

~~D 

150  RPM 

200  RPM 

0.125 

X 

X 

0.250 

X 

X 

0.500 

X 

X 

0.600 

X* 

1.000 

X 

X 

Propeller  3S36.  This  propeller  was 

built  and  tested. 

propellers  using  Lerbs*  induction  factors  have 
been  programmed  on  the  UNIVAC  high-speed 
digital  computer  at  the  Model  Basin. 10  The 
radial  distribution  of  wake  fraction  that  was 
used  in  the  propeller  calculations  is  presented 
in  Appendix  D  (see  Figure  20),  but  the  nominal 

U-^o) 

wake  was  corrected  by  the  -  and  re- 

d-^0.7) 

plotted  to  a  base  of  R/RQ  in  Figure  5. 

Curves  of  Ae/A0  versus  LQ/D  Are  shown 
in  Figure  23  for  N  =  150  and  200.  A  final  radial 
distribution  of  pitch  is  plotted  in  Figure  24 
(Appendix  E)  for  each  propeller  at  its  proper 
longitudinal  position  L0/D .  These  final  pitch 
curves  result  from  applying  appropriate  correc¬ 
tions  based  on  methods  presently  employed  at 


the  Model  Basin  for  friction  and  lifting-line  to 
lifting-surface  effects.8 
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Figure  5  —  Curves  of  Estimated  Effective  Velocity  Ratio  versus  R/R0 

PROPELLER  PERFORMANCE  AS  RELATED  TO  SPACING 


A  computed  curve  of  propeller  open-water  efficiency  versus  L0/D  is  given  in  Figure  6 
for  150  rpm  and  200  rpm.  The  propeller  efficiency  was  obtained  from  a  numerical  solution8,10 
of  Equation  [1];  see  Analysis.  Of  importance  is  the  variation  of  the  propulsive  coefficient 
T)p  with  L0/D ,  Figure  7  gives  an  estimated  curve  of  rjp  versus  L0/D  for  optimum  propellers 
where  the  components  in  the  relation  Tjp  »  ep  •  er  •  e A  are  obtained  in  the  following  manner: 
Propeller  efficiency  ( ep  .  er)  behind  the  body  is  approximated  by  Equation  [1]  using  local 


1  -  t 


wake  fractions  from  Figure  5.  Hull  efficiency  eh  =  - is  estimated  by  using  the  thrust 

1  -  H^o 

deduction  data  in  Figure  3  and  by  assuming  the  effective  wake  fraction  as  being  equal  to 
the  wake  taken  at  R/R0  *  0.7  from  Figure  5. 
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Figure  6  -  Propeller  Efficiencies  Calculated  According  to  Equation  ]la] 
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Figure  7  -  Propulsive  Coefficient  versus  L0/D 
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It  seems  logical  to  conclude  from  Figure  7  that  variations,  within  the  range  of  investi¬ 
gation,  in  the  longitudinal  position  of  the  propeller  do  not  critically  affect  propulsive  effici¬ 
ency.  In  fact,  it  appears  that  there  is  a  broad  flat  optimum  which  allows  considerable  freedom 
of  choice  in  locating  a  stern  propeller,  but  a  modest  increase  in  seems  likely  if  a  propeller 
is  positioned  somewhat  aft  of  the  reference  location. 


EXPERIMENTS 

To  experimentally  verify  the  calculated  performance  of  the  optimum  propellers, 
Propeller  3836  was  tested  in  open  water  and  behind  the  body  at  L0/D  «  0.6.  Appendix  B 
contains' a  description  of  the  salient  features  of  the  instrumentation  and  test  procedure. 
Resistance  tests  were  conducted  on  the  model  with  the  sting  extended  (see  Figure  8)  and 
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Figure  8  -  Model  4620  Equipped  for  Propulsion  Using  Propeller  3836 

with  the  sting  retracted.  No  measurable  change  in  the  residual-resistance  coefficient  Cr  was 
detected.  However,  it  is  believed  that  the  total  drag  would  change  in  proportion  to  the  wetted 
area  added  by  the  sting,  but  this  is  negligibly  small,  and,  therefore,  the  propulsive  coefficient 
can  be  used  as  a  powering  criterion  instead  of  shaft  horsepower. 

Open-water  characteristic  curves  for  Propeller  3836  are  shown  in  Figure  25  of  Appendix 
E.  For  L0/D  *  0.6,  the  corresponding  design  thrust  and  speed  coefficients,  Kt  and  J>  are 
0.192  and  1,292,  respectively.  Figure  25  shows  that  Propeller  3836  is  operating  at  practically 
maximum  efficiency  for  J  *  1.292.  The  experimental  propeller  efficiency  is  shown  in  Figure  6 
for  Propeller  383t>.  Based  on  these  data,  the  agreement  between  the  predicted  propeller  per¬ 
formance  and  the  experimental  result  is  good.  The  only  additional  determination  made  from 
the  open-water  data  is  the  effective  wake  fraction  as  derived  from  the  usual  Froude  synthesis. 

Turning  now  to  the  propulsion  problem,  first  compare  the  hull  efficiency  quantities 
(1  - 1 )  and  (1-W)  which  were  obtained  by  experiment  and  by  computation.  The  results  of  the 
submerged  propulsion  test  are  given  in  Figure  9  as  curves  of  the  various  propulsion  coefficients 
versus  the  total  resistance  coefficient.  Figure  9  is  entered  at  the  indicated  design  J  with  the 
following  result:  (1  -t)  =  1.00  and  (1  -W  )=  0.85.  Using  the  faired  values,  (1  -t)  *  0.98  from 
Figure  3  and  (1-IV0  7)  *  0-84  fr°m  Figure  5,  it  is  found  that  the  computed  and  experimental  hull 
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Figure  9  -  Submerged  Propulsion  Test  with  Propeller  3836 

efficiency  quantities  are  in  close  agreement.  This  close  agreement  between  the  computed  and 
experimental  results  for  the  hull  efficiency  is  considered  an  important  consequence  of  this 
study.  The  experimental  curve  of  rjp  in  Figure  9  was  determined  directly  from  measured  values 
of  model  resistance  and  propeller  power.  Although  the  er  curve  of  Figure  9  was  calculated  as 
a  derived  quantity,  the  experimental  value  eT  =  0.94  at  design  J  compares  well  to  a  computed 
value  er  *  0.95  from  Equation  [1]  and  [la].  The  experimental  value  of  r}p  is  about  0.91  at 
design  J .  This  value  is  plotted  in  Figure  7  at  L0/D  *  0.6,  where  it  is  seen  to  lie  on  the  com¬ 
puted  curve. 


CONCLUDING  REMARKS 

The  effect  on  the  propulsion  quantities,  thrust  deduction,  wake,  propeller  efficiency, 
and  propulsive  coefficient  of  variations  in  body-propeller  axial  clearance  at  deep  submergence 
on  a  body  of  revolution  Series  58  form  was  determined.  Some  principal  findings  which  have 


resulted  from  this  study  are:  (1)  The  pressure  field  due  to  propeller  action  is  significant  for 
approximately  the  last  20  percent  of  the  body  length,  and  the  thrust-deduction  coefficient  is 
essentially  zero  for  body-propeller  clearances  (measured  from  a  reference  station  at  0.966  of 
the  body  length)  greater  than  the  maximum  diameter  of  the  body.  (2)  A  rather  flat  curve  of 
propulsive^coefficient  versus  body- propeller  clearance  was  calculated,  but  a  modest  increase 
of  several  points  in  propulsive  coefficient  might  be  possible  if  the  propeller  is  located  at 
about  0.6  diameters  behind  the  reference  station.  For  the  condition  investigated,  the  experi¬ 
mental  results  are  in  good  agreement  with  the  computed  values.  For  example,  at  a  body- 
propeller  clearance  of  0.6  maximum  body  diameters  comparative  results  were  obtained  as 
follows: 


Quantity 

Computed 

Experiment 

Thrust  deduction  factor 

0.98 

1.00 

Wake  factor 

0.84 

0.85 

Propeller  open-water  efficiency 

0.82 

0.82 

Relative  rotative  efficiency 

0.95 

0.94 

Propulsive  coefficient 

0.91 

0.91 

Although  the  foregoing  discussion  has  dealt  with  propeller  location  in  connection  with  power¬ 
ing  aspects  of  a  bare  hull,  the  need  for  collateral  information  on  vibration,  cavitation,  and 
noise  is  emphasized.  For  example,  in  a  circumferentially  nonuniform  wake,  the  wake  peaks 
are  diminished  by  increasing  the  axial  clearance  between  an  appendage  configuration  and  the 
propeller.  Two  possible  outcomes  of  increased  clearance  are  (1)  improved  propeller  cavitation 
and  noise  characteristics,  and  (2)  lower  fluctuating  thrust  amplitudes.  An  additional  improve¬ 
ment  which  might  occur  as  a  result  of  increased  hull-propeller  clearance  is  a  reduction  of 
surface  forces  due  to  hydrodynamic  pressure  fluctuations.11  This  is  an  important  problem  in 
connection  with  vibration  and  noise. 
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APPENDIX  A 

GEOMETRICAL  DATA  AND  PHOTOGRAPH  OF  MODEL  4620 
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TABLE  2 


Offsets  and  Particulars  for  Model  4620 


X 

y 

X 

y 

in. 

in. 

in. 

in. 

0.0 

0.0 

93.6 

11.82 

3.6 

3.5 

97.2 

11.66 

7.2 

4.977 

100.8 

11.49 

10.8 

6.108 

104.4 

11.29 

14.4 

7.047 

108.0 

11.07 

18.0 

7.850 

111.6 

10.83 

21.6 

8.549 

115.2 

10.56 

25.2 

9.160 

118.8 

10.27 

28.8 

9.697 

122.4 

9.954 

32.4 

10.17 

126.0 

9.613 

36.0 

10.58 

129.6 

9.243 

39.6 

10.93 

133.2 

8.843 

43.2 

11.24 

136.8 

8.411 

46.8 

11.50 

140.4 

7.945 

50.4 

11.71 

144.0 

7.447 

54.0 

11.89 

147.6 

6.910 

57.6 

12.03 

151.2 

6.334 

61.2 

12.13 

154.8 

5.715 

64.8 

12.21 

158.4 

5.053 

68.4 

12.25 

162.0 

4.344 

72.0 

12.27 

165.6 

3.584 

75.6 

12.25 

169.2 

2.774 

79.2 

12.21 

172.8 

1.908 

82.8 

12.15 

176.4 

0.984 

86.4 

12.06 

180.0 

0.0 

90.0 

11.97 

Model  Particulars: 

Length,  ft 

*  15.0 

Diameter,  ft 

=  2.044  (24.53  in.) 

Nose  radius,  ft 

0.1392  (1.670  in.) 

Tail  radius,  ft 

=  0.0 

Vetted  surface,  ft^ 

=  70.55 

Volume,  ft^ 

=  29.53 

LCB,  ft 

6.6840 
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APPENDIX  B 

INSTRUMENTATION  AND  TEST  PROCEDURE 
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WAKE  SURVEYS 


Because  of  the  flow  symmetry  without  appendages  a  3-dimensional  flow  survey  was 
felt  to  be  unnecessary  and,  consequently,  the  arrangement  of  cylindrical  total  head  and  static 
head  tubes  shown  in  Figures  12  and  13  was  used.  Six  transverse  surveys  were  made;  the 
longitudinal  position  of  each  plane  of  measurement  was  as  follows:  L0/D  *  0,  L0/D  *  1/8, 
Lq/D  *  1/4,  Lq/D  ■  1/2,  Lq/D  =  1.0,  and  L0/D  ■  1.5.  Except  for  L0/D  *  0,  the  surveys 
extend  radially  to  R/RB  -  0.53,  where  RB  =  1.022  feet.  At  Station  L0/D  *  0,  a  spacer  block 
inserted  at  the  mounting  head  increased  R/RB  to  0.571.  This  adjustment  was  made  to  provide 
a  minimum  clearance  of  four  tube  diameters  between  the  hull  surface  and  the  innermost  pitot 
tube.  The  entire  pressure-measuring  system  consisted  of  the  pitot  rakes  as  the  sensing 
device  and  straight  tube  manometers  for  the  measuring  element.  The  use  of  water  as  the 
metering  fluid  allowed  an  accuracy  of  about  0.15  percent  of  stagnation  pressure  to  be  obtained 
at  the  model  test  speed  of  8  knots. 

The  pitot  rakes  were  connected  to  straight-tube  manometer  boards  mounted  on  the 
towing  carriage.  A  large  capacity  vacuum  tank  was  connected  in  parallel  to  a  manifold 
located  at  the  top  of  each  board,  and  a  suitable  vacuum  was  applied  to  each  manifold  to 
establish  a  board  level.  A  reference  tank  partially  filled  with  water  and  open  to  atmospheric 
pressure  was  connected  to  a  reference  tube  on  each  manometer  board.  A  horizontal  datum 
line  was  formed  by  the  zero  reading  (model  at  rest)  of  each  tube.  For  each  board,  the  dif¬ 
ference  between  the  reference  tube  zero  reading  and  the  reference  reading  at  the  equilibrium 


Figure  12  —  Pitot  Rakes  as  Mounted  on  Model  4620 


run  was  applied  as  a  correction  to  tbe  datum  line,  to  account  for  the  effect  on  the  aero  readings 
of  vacuum  changes  in  tbe  system.  The  equilibrium  run  was  determined  by  observing  the  maxi¬ 
mum  or  minimum  heads  attained  from  a  series  of  constant  speed  runs.  Normally,  about  five 
runs  were  necessary  to  reach  equilibrium.  Opening  and  closing  the  valves  of  each  tube  was 
accomplished  by  actuating  a  solenoid  valve  immediately  after  and  before  the  acceleration  and 
deceleration  parts  of  a  run,  respectively. 

RESISTANCE  AND  POWERING  EXPERIMENTS 

Resistance  tests  were  conducted  by  towing  the  model  with  the  single  strut  and  internal 
resistance  dynamometer  arrangement  shown  by  Figure  14.  With  the  model  still  rigged  on  the 
single  strut,  propulsion  tests  using  a  submersible  motor  and  propulsion  dynamometer  were 
made  immediately  following  the  resistance  test.  This  procedure  ensured  the  same  orientation 
and  condition  of  the  model  for  the  two  tests.  Both  resistance  and  propulsion  tests  were  con¬ 
ducted  with  a' turbulence-stimulating  sand  strip  placed  around  the  girth  of  a  body  section  at 
1/20  of  the  body  length.  When  necessary,  appropriate  corrections  were  applied  to  the  drag 
data  for  towing-strut  interference  effects  and  for  the  parasitic  drag  of  the  sand  strip.  The 
submerged  propulsion  tests  were  conducted  at  constant  model  speed.  Test  speeds  of  6  and 
8  knots  were  selected  in  order  to  cover  the  desired  range  of  propeller  loading  within  the  limits 
of  dynamometer  capacity.  The  propeller  was  driven  by  a  5-hp  motor  and  a  transmission-type 
dynamometer,  with  a  capacity  of  100  lb  of  thrust  and  100  in-lb  of  torque,  placed  between  the 
propeller  and  the  driving  motor. 

In  essence,  the  various  propulsion  quantities  are  obtained  as  functions  of  Cg  much  in 
the  same  manner  as  open-water  propeller  characteristics  are  obtained  as  functions  of  J,  In 
both  cases  the  data  are  expressed  as  functions  of  load. 
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Figure  14  -  Schematic  Diagram  of  Arrangement  of  Model-Towing  Apparatus 
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APPENDIX  C 

DATA  FOR  USE  IN  COMPUTING  THRUST  DEDUCTION  FOR  MODEL  4620 
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THRUST  DEDUCTION  COMPUTATIONS 


The  nondimensional  disturbing  velocities  v/va  and  v^/u0,  where  utf  is  the  propeller- 
induced  velocity  at  infinity  due  to  a  cylindrical  vortex  sheet*  of  ring  vortices,  are  given  as 
functions  of  the  spherical  coordinates  r/R  and  6  in  Reference  1.  These  velocities  are  used 
to  obtain  the  velocity  ratio  V/V .  The  procedure  for  calculating  V/V  is  outlined  briefly.  A 
radial  distribution  of  thrust  or  circulation  required  for  the  calculation  of  thrust  deduction  for 
a  wake-adapted  propeller  can  be  approximated  by  an  equivalent  rectangular  distribution.  Con¬ 
sistent  with  this  approximation  and  the  characteristics  of  the  mathematical  model  of  the  pro¬ 
peller  (i.e.,  infinite  number  of  blades  and  very  small  pitch),  it  has  been  shown  in  Reference  1 


that  the  rate  of  change  in  the  2-direction  of  the  elementary  circulation  —  about  a  vortex 

dX 

filament  is  equal  to  the  induced  velocity  ua  far  behind  the  propeller  disk*  In  other  words,  the 
circulation  per  unit  length  at  a  sufficiently  large  distance  behind  the  propeller  disk  is  u  . 
Numerically,  write 


From  momentum  theory,1  velocity  ua  is  related  to  propeller  thrust  coefficient  CT  and  the  pro¬ 
peller  speed  of  advance  Va  by  the  quadratic  equation 


which  when  solved  for  ua/VQ  becomes 


and 


uo  _ 

_  ..-1+>/l+Cr 

<U 


[3] 


[3a] 


where  V  is  the  body  speed  and  W0  is  the  effective  wake  fraction* 

It  has  been  shown1  that  a  representative  radial  thrust  distribution  is  given  by  a  pair  of 

er 

cylindrical  vortex  sheets,  one  at  R/R0  »  0.23  with  a  circulation  distribution - and  the 

„  dX 

dr 

other  at  R/R0  -  0.85  with  strength  —  .  The  effective  value  of  CT  is  based  on  the  reduced 


♦This  cylindrical  vortex  sheet  is  composed  of  a  succession  of  closely  spaced  circular  vortex  filaments  which 
extend  to  infinity  behind  the  propeller. 


area;  the  ratio  r/R,  on  the  reduced  radius.  The  procedure  for  calculating  the  velocity  ratio 
U/V  is  now  complete  with  the  finding  of  ua/V ;  it  is  possible  to  express  all  velocities  as 
fractions  of  the  speed  V.  The  pressure  coefficients  are  obtained  from  the  Bernoulli  equation 


where  q  is  the  free-stream  stagnation  pressure.  The  thrust»deduction  coefficient  t  is  a 
function  of  the  integral  of  the  difference  A p/q  between  the  towed  and  propelled  pressure 
coefficients  calculated  for  the  afterbody.  When  the  pressure  coefficient  is  invariant  around 
the  circumference  of  a  section  (bare-hull  case),  it  can  be  shown  that 


where 


and 


T  -  Rt  AS,  n02q 
T  "  T  "  T 


A Rt  -  (ir/?2?)ACD  -  *D2q 


b  —  f\x)dz 
? 

A  —  f\x)dx 
? 


dy 


1  dA 
~vD2  dx 


[4] 


The  fuQction  f\x)  is  the  slope  of  the  nondimensional  sectional  area  curve. 
A  Series  58  form12  (TMB  Model  4620) 


was  selected  for  the  computable  example  and 
for  the  experiments.  Model  particulars  are 
given  in  Appendix  A.  A  reference  station 
L0  -  0  (see  Figure  1  in  Analysis),  is  located 
at  0.966  of  the  body  length.  The  function 

2 y-— ^  has  been  computed  for  the  last  20  p®*- 
dx 

cent  of  the  body  length,  and  its  values  are 

tabulated  in  Table  8. 

Given  the  geometry  of  the  body,  the 
thrust  deduction  coefficient  t  is  computed  by 
using  Equation  [$a]f  the  propeller  nondimen- 
sional  disturbing  velocities  v/va  *nd 
and  Equation  [4].  The  propeller  designs  for 
the  various  longitudinal  positions  L^/D  are 


TABLE  3 

Slope  of  Nondimensional  Sectional  Area 
Curve  for  Series  58  Form,  Model  4620 


X 

-§*|3 

3* 

CM 

0.80 

-0.64078 

0.82 

-0.63870 

0.84 

-0.62828 

0.86 

-0.60792 

0.88 

-0.57574 

0.90 

-0.52966 

0.92 

-0.46741 

0.94 

-0.38632 

0.96 

-0.28350 

0.98 

-0.15590 

1.00  • 

0.00000 

for  a  constant  thrust  T  of  82,680  lb.  An  effective  CT  for  Equation  [8]  is  calculated  for  each 
longitudinal  position  by  obtaining  V  from  the  typical  resistance  curve  depicted  in  Figure  15. 
In  Figure  16  the  curve  was  entered  on  the  ordinate  Rt  *  82,680  (1-/),  where  (1  -t)  was  esti¬ 
mated  from  data  of  Breslin2  and  available  experimental  data  at  L0/D  «  1/8;  see  data  point  in 
Figure  8..  The  results  of  this  procedure  are  summarized  as  follows: 


(1-0 

d-%)* 

V 

knots 

C 

0.125 

0.920 

0.805 

30.0 

0.5118 

0.25 

0.950 

0.817 

30.42 

0.4834 

0.50 

0.978 

0.837 

30.81 

0.4488 

1.0 

1.00 

0.865 

31.09 

0.4127 

*  Appropriate  values  of  ( 1  —  were  taken  from  Figure  5. 


Since  the  methods  and  computer  programs  available  for  determining  the  potential  flow  at  the 
body  surface  for  the  body  alone  are  considered  elsewhere,13  they  will  not  be  discussed  here. 
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Figure  15  -  Typical  Resistance  Curve 
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TABLE  5 
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14.  Column  13  —  Column  10  and  Equation 


s  a  b> 


1=  hr  b. 


=  hr  b> 


00  h>  b- 


?  I  3° 


Cvi  ?>> 


Hull  Propeller  Total 


0.80  0.3036  7.799  8.182  10.527  0.002256  0.000  0.000368  0.000  1.00000  0.000368  1.0004  -0.0008  0.0008 

0.88  0.2060  7.205  10.780  7.7281  0.004186  0.000683  0.9716  0.000682  0.9723  0.0546  0.0014 

0.92  0.1461  6.241  12.295  6.3228  0.006254  0.001020  0.9487  0.00100  0.9497  0.0981  0.0019 

0.96  0.0778  4.269  13.942  4.9179  0.01034  0.001687  0.8944  0.00166  0.8961  0.1970  0.0030 

0.98  0.0401  2.564  14.835  4.2175  0.01406  0.002294  0.8031  0.00224  0.8053  0.3515  0.0035 
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APPENDIX  D 
WAKE  DATA 
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PITOT  RAKE  CALIBRATION 


The  size  of  the  tubes,  the  radial  spacing  between  tubes,  and  the  form  of  the  fairing  are 
in  accordance  with  design  practice;  see  Appendix  B.  In  contrast  to  spherical-headed  pitot 
tubes  which  are  sensitive  to  turbulence  it  appears  that,  within  the  range  of  the  test,  Reynolds 
number  has  a  rather  small  effect  on  the  readings  obtained  from  a  pitot  rake  of  the  type  used  in 
the  present  study.  Calibration  data  were  obtained  at  two  carriage  speeds,  8  and  6  knots.  The 
results  of  the  calibration  tests  are  shown  in  Figures  16  and  17  where  the  net  impact  and  the 
net  static  heads  are  plotted  versus  R  for  each  speed.  In  the  relation 

V  =  Va  = 

the  coefficient  C  was  calculated  to  be  1.01  for  the  6-knot  test  and  1.006  for  the  8-knot  test. 

An  average  value  of  C  «  1.008  was  used  to  reduce  all  subsequent  wake  data  behind  the  body. 

DISCUSSION  OF  RESULTS  OF  WAKE  SURVEYS 

Free-turbulence  shear-flow  theory  and  various  similarity  hypotheses  are  not  valid  in 
regions  of  considerable  pressure  gradients.  The  wake  traverses  obtained  at  the  six  stations 

(Lq/D  =  0  to  Lq/D  =  1.5)  are  essential  to  the  solution  of  the  problems  of  this  study  and  will 
undoubtedly  be  of  additional  use  in  other  problems.  The  original  test  data  for  each  wake  sur¬ 
vey  showed  slight  scatter  when  plotted  as  curves  of  total  head  and  static  head  versus  the 
radius  R.  Within  the  limits  of  test  data,  the  variation  of  static  head  with  R  was  linear.  For 
purposes  of  analysis  and  discussion,  the  original  data  were  crossfaired  and  replotted  in  non- 
dimensional  form  in  Figures  18,  19,  and  20.  The  most  interesting  features  of  the  curves  shown 
in  Figure  18  are:  (1)  The  asymptotic  nature  of  the  wake  curves  beginning  at  about  L0  *  D>  and 
(2)  the  relatively  small  change  with  LQ/D  at  R/Rg  >  0.5. 

In  Figure  19,  the  variation  of  the  pressure  coefficient  p/q  with  LQ/D  is  shown  with 
R/Rg  as  a  parameter.  It  is  clear  from  these  curves  that  the  static  gradient  is  in  the  direction 
of  increasing  radius;  i.e.,  at  each  position  LQ/£>,  the  velocity  is  increasing  with  radius.  This, 
of  course,  is  the  expected  behavior;  however,  the  real  significance  of  the  curves  in  Figure  19 
is  the  manner  in  which  the  p/q  curves  converge  at  higher  values  of  L0/D .  At  LQ  *  Dy  there 
is  very  little  spread  with  the  parameter  R/Rg  and  at  L0/D  *  1.5,  the  static  pressure  is  invariant 
radially  within  the  range  of  measurements.  One  advantage  of  conducting  a  wake  survey  by 
taking  separate  measurements  of  total  head  and  static  head  is  the  ability  to  examine  the  static 
gradient  and  to  calculate  the  potential  component  of  the  wake. 

CALCULATIONS  FOR  POTENTIAL  WAKE 

The  streamline  wake  can  be  calculated  from  pitot-static  measurements  based  on  the  fol¬ 
lowing  considerations;  (1)  The  fluid  is  frictionless,  no  energy  loss  occurs,  and  (2)  the  model 
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Figure  16  -  Pitot  Hake  Calibration  Data  at  6  Knots 
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Figure  17  -  Pitot  Rake  Calibration  Data  at  8  Knots 
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Figure  19  —  Distribution 
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Figure  20  —  Distribution  of  (1  -W)  as  a  Function  of  R/RB 
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is  steadily  advancing  in  undisturbed  water.  When  -  0,  the  Bernoulli  constant  as  measured 
by  a  pitot  tube  at  some  distance  ahead  of  the  model  is  q  so  that 

%P(V-Vp)2+p-q  [6] 

where  Vp  -  VWp,  the  potential  wake  velocity.  Finally,  in  terms  of  the  pressure  coefficient, 
the  potential  velocity  ratio  becomes 

[6] 

v  9 

Values  of  (1  -Wp)  were  calculated  from  Equation  [6]  and  are  plotted  in  Figure  21a.  A  slight 
spread  in  (1  )  with  the  parameter  R/RB  is  observed  in  the  region  of  L0/D  <  0.25.  Maxi¬ 

mum  change  in  Wp,  within  the  range  of  the  test  conditions,  is  in  the  order  of  five  points.  This 
corresponds  to  the  small  displacement  flow  associated  with  well  streamlined  bodies.  The 
single  curve  of  (l-Wp )  (see  Figure  21a)  at  L0/D  -  1.5  will  decelerate  and  approach  unity  at 
some  greater  distance  downstream;  the  total  velocity  ratio  (1— IV)  will  also  become  unity,  but 
will  not  reach  free-stream  velocity  until  a  much  greater  distance  downstream  than  the  ratio 
(1-1V  ).  In  fact,  the  concept  that  the  potential  flow  field  decays  at  a  greater  rate  than  does 
the  frictional  flow  field  is  a  basis  for  s  peculating  that,  from  the  propulsion  viewpoint,  an 
optimum  propeller  location  exists. 

Figure  21b  shows  curves  of  potential  velocity  ratio  for  the  case  of  pure  streamline  flow. 
This  figure  also  shows  that  the  velocities  at  equal  values  of  L0/D  are  several  points  lower  than 
those  indicated  by  the  comparable  curves  of  Figure  21a.  Although  the  curves  of  Figure  21b  are 
mostly  pt  academic  interest,  several  points  concerning  the  theoretical  calculation  of  (1-Wp) 
deserve  discussion.  Calculation  of  a  “so-called**  potential  velocity  ratio  from  experimental 
measurements  is  a  scheme  for  resolving  the  total  velocity  of  a  viscous  flow  into  component 
parts.  The  procedure  followed  in  the  derivation  of  the  curves  presented  in  Figure  21a  is 
really  an  attempt  to  evaluate  the  velocity  due  only  to  normal  stress.  Of  course,  the  wake 
behind  the  hull  is  a  free  turbulent  shear  flow.  Referring  now  to  Figure  21b,  the  results  shown 
are  for  a  closed  body  boundary  condition.  An  exact  general  method  of  solving  the  Neumann 
problem  has  been  developed  for  axisymmetric  flow.  The  solution  is  programmed  on  a  high¬ 
speed  computer  and  solutions  can  be  obtained  for  both  closed-  and  open-body  conditions.  In 
an  effort  to  approximate  the  mathematical  representation  of  the  real  wake  behind  the  hull,  a 
new  body  shape  with  a  tail  was  derived  by  adding  the  displacement  thickness  of  the  boundary 
layer  to  the  original  body.  Agreement  with  the  experimental  results  was  not  improved  by  this 
process,  since  one  undesirable  feature  of  the  body  plus  tail  condition  is  the  exclusion  of 
points  close  to  the  centerline  behind  a  body.  Of  course,  the  greatest  velocity  change  is 

within  this  very  region. 
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Figure  23  -  Expanded-Area  Ratio  versus  LJD  for  Optimum  and  Nonoptimum  Propellers 


Finol  Pitch  Rotio 


Figure  24  -  Final  Radial  Distribution  of  Propeller  Pitch  Ratio  for  Each  Condition 
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THRUST  COEFFICIENT  (Kt),  TORQUE  COEFFICIENT  (lOKq)  AND  EFFICIENCY  (•) 
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30  MARCH  1960 
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SPEED  COEFFICIENT  J 


Figure  25  —  Open-Water  Characteristic  Curves  for  Propeller  3836 


46 


REFERENCES 


1.  Korvin-Kroukovsky,  B.V.,  “Stern  Propeller  Interaction  with  a  Streamline  Body  of 
Revolution,”  International  Shipbuilding  Progress,  Vol.  3,  No.  17  (Jan  1956). 

2.  Breslin,  J.P.,  “A  Simplified  Theory  for  the  Thrust  Deduction  Force  on  a  Body  of 
Revolution,”  Proceedings  5th  Midwestern  Conference  on  Fluid  Mechanics  (1957). 

3.  Astrup,  N.C.,  “A  Set  of  Propellers  in  Open  Water  and  behind  a  Body  of  Revolution,” 
Parts  1  and  2,  Norwegian  Ship  Model  Experiment  Tank  Publication  Nos.  45  and  46  (1957). 

4.  Hickling,  R.,  “Propellers  in  the  Wake  of  an  Axisymmetric  Body,”  Institution  of  Naval 
Architects,  Quarterly  Transactions,  Vol.  99,  No.  4,  (Oct  1957). 

5.  Weinblum,  G.,  “The  Thrust  Deduction,”  American  Society  of  Naval  Engineers,  Vol.  63 
(1951). 

6.  Dickmann,  H.E.,  “The  Interaction  between  Propeller  and  Ship  with  Special 
Consideration  to  the  Influence  of  Waves,”  Jahrbuch  der  Schiffbautechnischen  Gesselschaft, 

40  (1939). 

7.  Schoenherr,  K.E.,  and  Aquino,  A. A.,  “Interaction  between  Propeller  and  Hull,”  U.S. 
Experimental  Model  Basin  Report  470  (Mar  1940). 

8.  Eckhardt,  M.K.  and  Morgan,  W.B.,  “A  Propeller  Design  Method,”  Transactions  Society 
of  Naval  Architects  and  Marine  Engineers,  Vol.  63  (1955). 

9.  Gawn,  R.W.L.,  “Effect  of  Pitch  and  Blade  Width  on  Propeller  Performance,” 
Transactions  Institution  of  Naval  Architects,  Vol.  95  (1953). 

10.  Hecker,  Richard,  “Manual  for  Preparing  and  Interpreting  Data  of  Propeller  Problems 
Which  Are  Programmed  for  the  High-Speed  Computers  at  the  David  Taylor  Model  Basin,” 

David  Taylor  Model  Basin  Report  1244  (Aug  1959). 

11.  Tachmindji,  A.J.  and  Dickerson,  M.C.,  “The  Measurement  of  Oscillating  Pressure  in 
the  Vicinity  of  Propellers,”  David  Taylor  Model  Basin  Report  1130  (Apr  1957). 

12.  Landweber,  L.  and  Gertler,  M.,  “Mathematical  Formulation  of  Bodies  of  Revolution,” 
David  Taylor  Model  Basin  Report  719  (Sep  1950). 

13.  Smith,  A.M.O.  and  Pierce,  J.,  “Exact  Solution  of  the  Neumann  Problem,  Calculation 

of  Non-Circulatory  Plane  and  Axially  Symmetric  Flows  about  or  within  Arbitrary  Boundaries,” 
Douglas  Aircraft  Company  Report  ES  26988  (Apr  1958). 


47 


INITIAL  DISTRIBUTION 


Copies 

9  CHBUSHIPS 

3  Tech  Lib  (Code  210L) 

1  Lab  Mgt  (Code  320) 

2  Prelim  Des  Br  (Code  420) 

1  Mach,  Sci,  &  Res  Sec  (Code  436) 

1  Sub  Br  (Code  525) 

1  Prop,  Shaft,  &  Bear  Br  (Code  644) 

S  CHBUWEPS 

2  Tech  Lib  (DL  1-3) 

1  Fluid  Mech  &  Flight  Dyn  Sec  (RRRE-4) 
1  Supporting  Res  Br  (RUTO-3) 

1  RUTO-32 

Attn:  Mr.  H.  Eggers 

3  CHONR 

1  Fluid  Dyn  Br  (Code  438) 

1  Struc  Mech  Br  (Code  439) 

1  Undersea  Programs  (Code  466) 

1  CDR, USNOL 

1  DIR,  USNRL 

1  CO  &  DIR,  USNEES 

1  CDR,  U SNOTS,  Pasadena 

1  DTMB  High-Speed  Phenomena  Div 
Langley  Field 

1  SUPT  USNAVP6SC0L 

1  NAVSHIPYD  PTSMH 

1  NAVSHIPYD  MARE 

1  0  in  C,  USN  Sub  School 

Attn:  Sub  Dept 

1  0  in  C,  PGSCOL,  Webb 

1  DIR,  Aero  Res,  NASA 

10  CDR,  ASTIA 

1  Dir  Def  R  &  E 

1  Hydro  Lab,  CIT 

2  Dept  NAME,  MIT 

1  DIR,  St.  Anthony  Falls  Hydrau  Lab 

1  DIR,  ORL 

Attn:  Dr.  G.F.  Wisliconus 

1  DIR,  DL,  SIT 

2  Iowa  Inst  of  Hydrau  Res 

State  Univ  of  Iowa 

1  Dr.  Hunter  Rouse 
1  Dr.  L.  Landweber 


Copies 

2  Dept  NAME,  Univ  of  Michigan 

Attn:  Prof  R.B.  Couch 

1  Univ  of  Notre  Dame 

Attn:  Dr.  A.G.  Strandhagen 

3  School  of  Engin 
Univ  of  California 

1  Prof  H.A.  Sdude 
1  Dept  of  Naval  Arch 
1  Engr  Libr 

1  Oceanics,  Inc. 

Attn:  Dr.  Paul  Kaplan,  Pres 

1  TRG 

1  Hydronautics,  Inc. 

1  Ingalls  Shipbldg  Corp. 

1  New  York  Shipbldg  Corp. 

2  Gen  Dyn  Corp.,  EB  Div. 

1  Bethlehem  Steel  Co.,  Shipbldg  Div.,  Quincy 

Attn:  Mr.  Hollinshead  DeLuce 

2  Gibbs  &  Cox,  Inc. 

1  NNSB  &  DO  CO.,  Engin  Tech  Div 
Attn:  Mr.  John  Kane 

1  SNAME 

Attn:  Capt  W.N.  Landers 

1  Douglas  Aircraft  Co.,  Inc. 

Aircraft  Div.,  Long  Beach 
Attn:  Mr.  John  L.  Hess 

1  Vidya,  Inc. 


49 


